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Abstract 

In this paper we consider the problem of determining an unknown 
pair A, /i of piecewise constant Lame parameters inside a three dimen- 
sional body from the Dirichlet to Neumann map. We prove uniqueness 
and Lipschitz continuous dependence of A and \x from the Dirichlet to 
Neumann map. 

1 Introduction 

A relevant inverse problem arising in nondestructive testing of materials is the 
one of determining, within an isotropic, linearly elastic three dimensional body 
the elastic properties of the body from measurements of traction and dis- 
placement taken on the exterior boundary of the domain Q. 
This leads mathematically to the formulation of the following boundary value 
problem for the system of linearized elasticity 

div(CVu) = in ft C M. 3 , (1) 
u = ip on <9£1, 

where f2 is an open and bounded domain, V« denotes the strain tensor Vm := 
| (Vu + (Vu) T j, ip G H 1 ^ 2 (dfl) is the boundary displacement field, and C G 
L°°(f2) denotes the isotropic elasticity tensor with Lame coefficients A, /i: 

C = A/3 ® I3 + 2{jl sym , a.e. in 0, 

where -Z3 is 3 x 3 identity matrix and Isym is the fourth order tensor such that 
^SymA = A, The strong convexity condition is assumed 

A 1 > > 0, 2fi + 3A > p > a.e. in fl. 
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Under the above assumptions problem ([T]) has a unique weak solution u G H 1 (f2) 
and the Dirichlet-to- Neumann linear map (DN map), Ac, is well defined 

A c : ip e H 1/2 (ft) -> (CVn)^| a n e H^^idVi) 

where v is the exterior unit normal to d£l. 

The inverse problem consists in determining C, i.e. A and fi, from knowledge of 
the DN map A c . 

This problem is closely related to the conductivity inverse problem arising in the 
modelling of EIT (Electrical Impedence Tomography). For the mathematical 
treatment of this problem we refer to the fundamental papers [All] . |AP] . |Naj 
and [STJ) . 

Unfortunately, the mathematical approach used to investigate the conductivity 
inverse problem fails partly when dealing in the elasticity framework. This is 
due to the fact that we have to deal with an elliptic system instead that with a 
scalar equation and we have to recover two parameters A and \i instead of the 
sole conductivity parameter. As a consequence of these difficulties only partial 
results to this inverse problem are known and mainly concern the uniqueness 
issue. More precisely, the study of the problem was initiated in the 90's by 
Ikeata in [Ik] who considered a linearized version of it. In two dimensions Aka- 
matsu, Nakamura and Steinberg in )ANS] and for dimension n > 3 Nakamura 
and Uhlmann in |NU] showed that one can determine uniquely and in a stable 
way C°°(f2) Lame parameters A and fi and their derivatives on the boundary of 
a smooth domain fi from the DN map. Local uniqueness has been proved in di- 
mensions two by Nakamura and Uhlmann in |NU1] for C°°(f2) Lame parameters 
assuming that they are both close to positive constants. In three dimensions 
and higher Nakamura and Uhlmann in |NU2j and Eskin and Ralston in |ERj 
proved local uniqueness for smooth Lame parameters whenever fj, is close to a 
constant. 

To our knowledge no result concerning stability estimates is known. Based on 
the results obtained by Alessandrini in [All) who proved, logarithmic stability 
estimates for the conductivity inverse problem in the case of smooth conduc- 
tivities and the example of Mandache in [M] who proved the optimality of this 
estimate also for the inverse elasticity problem in the case of smooth Lame pa- 
rameters logarithmic stability estimates or even worse ones are expected. 
These considerations lead in recent years to look for different a priori assump- 
tions on the unknown parameters which take into account the applied context 
from which the problem arises and give rise to better stability estimates ([V], 
[ABV] . [BF] . [BFV] . [ABF],(BdHQ]). An attempt in this direction has been 
done by Alessandrini and Vessella in }AVj with for the conductivity inverse 
problem for unknown conductivities depending only on a finite number of pa- 
rameters. In |AV] they proved Lipschitz continuous dependence from the DN 
map for conductivities that are constant on known subdomains, assuming ellip- 
ticity and C 1,a regularity at the interfaces joining contiguous domains. 
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In this paper we propose to consider L°° (fl) elasticity tensors of the form 

N 
3=1 

where the Dj' s, j = 1, • • • , N, are known disjoint Lipschitz domains represent- 
ing a partition of fl and Xj,fJ,j, ]=!,■■■ ,N, are unknown constants. 
We will prove that if C 1 and C 2 are of this form, assuming that the boundaries 
of Dj's and of Q contain flat portions, we have 

IIC 1 - C 2 ||oo < C||A c i - ^c 2 \\c(HU2( d n),H-^ 2 (dn))- 

where the stability constant C appearing in the estimate depends on various 
parameters like ao , /?o , on the regularity bounds on 51 and on the Dj ' s and on 
their number N. In particular, the constructive character of the proof allows us 
to establish an estimate from above of the constant C. 

We want now to emphasize that several significant examples fit in our analy- 
sis. Polyhedral partitions of f2, appearing in any finite-element scheme used for 
effective reconstruction of the Lame parameters (see for example [BJK ; and a 
layered configuration of the sets Dj 's arising in the study of composite laminates 
|Mi] and in geophysical prospection, |BCj . 

Our approach is based on the use of the following key ingredients: existence of 
singular solutions and study of their behaviour close to the flat discontinuity 
interfaces of the Dj's, regularity estimates and quantitative estimates of unique 
continuation of solutions to system (p}. 

As already pointed out in |BF] a relevant difference with respect to the scalar 
case treated in |AVj is the issue of existence of singular solutions. In fact, in 
the case of strongly elliptic systems with L°° coefficients in dimension n > 3, 
existence of the fundamental solution and of the Green's function cannot be in- 
ferred without additional assumptions. In |HK] Hofmann and Kim prove their 
existence under the additional information that weak solutions of the system 
satisfy Dc Giorgi-Nash type local Holder estimates. It is clear that in the case 
of a polyhedral partition of fl solutions might not enjoy Holder regularity at 
edges. On the other hand, in order to obtain our result, it is enough to con- 
struct singular solutions and analyze their behavior in a Lipschitz subset JC at 
a given positive distance from edges. Nevertheless, while for the scalar (even 
for the complex valued treated in |BFj ) conductivity equation is fairly easy to 
get a solution of the equation close to a flat interface by using the fundamental 
solution for the Laplace equation and suitable reflection arguments, this seems 
not to be possible for the Lame system. In order to construct singular solutions, 
we make use of special fundamental solutions constructed by Rongved in [R] for 
isotropic biphase laminates. Furthermore, looking at solutions of the elasticity 
system in K, we can use the results of [CKVCj and of |LN| deriving regularity 
estimates for the solutions which allow us to obtain Holder estimates of unique 
continuation in K,. 

We would like to point out that our proof is based on the use of solutions having 
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boundary displacement fields supported in the flat portion Si of dft. Hence, 
our stability result also holds replacing the full DN map with the local DN map 
that we will define in Section 2. Moreover, we expect to derive similar stability 
estimates also the case of domains Dj's with C 1,Q portions of interfaces and this 
analysis will be object of a forthcoming publication. 

Finally, we would like to emphasize that Lipschitz stability estimates have be- 
come crucial also for the effective reconstruction of the unknown coefficients. In 
fact, recently, in |dHQSl| and |dHQS2| , de Hoop, Qiu and Scherzer have shown 
that Lipschitz type stability estimates imply local convergence of iterative re- 
construction algorithms and the radius of convergence of the iterates depends 
on the stability constant and hence an explicit determination of the dependence 
of such constant from the a-priori parameters, in particular from the partition 
number N, is crucial. 

The plan of the paper is the following: section [5] contains the description of 
the main result. In paragraph 12 . 1 1 we introduce the notation and main defini- 
tions. In paragraph 12.21 we state the main a priori assumptions and the main 
result. We also reformulate the inverse elasticity problem in terms of the non- 
linear forward map F acting on a finite-dimensional subset of R 2N and recall 
a general result (Proposition 12. 5[) that will let us show that F has a Lipschitz 
continuous inverse. 

Section [3] contains some auxiliary results. In particular, in paragraph [32] we 
collect some properties concerning the fundamental solution in biphase elastic 
isotropic materials introduced by Rongved in [R] and we prove existence of sin- 
gular solutions for our Lame system. In paragraph 13 . 31 we state some regularity 
results and estimates of unique continuation concerning solutions of piecewise 
constant Lame systems. 

In section|4]we give the proof of our main result by verifying that the forward 
map F corresponding to our inverse problem, satisfies all the assumptions of 
Proposition 12.51 thus concluding the proof. 

Finally, in the Appendix we recall some known results concerning solutions 
to Lame systems with constant coefficients and the proof of the estimate of 
unique continuation stated in section [3] 
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2 Main result 

2.1 Notation and main definitions 

For every x € R 3 we set x — (x' \xz) where x' € R 2 and X3 <E K. For every 
x £ R 3 , r and L positive real numbers we will denote by B r (x), B' r (x') and 
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Qr,h{x) the open ball in M 3 centered at x of radius r, the open ball in M 2 centered 
at x' of radius r and the cylinder B' r (x') x (x 3 — Lr,x 3 + Lr), respectively; 
in the sequel B r (0), B' r (0) and Q r ,i(0) will be denoted by B ri B' r and Q r ,L, 
respectively. We will also denote by R 3 + = {(x',x 3 ) e M 3 : x 3 > 0}, M'i = 
{(x',x 3 ) e M 3 : x 3 < 0}, B+ = B r n 1R 3 , and 5,7 = B r n R 3 . 
For any subset D of K 3 and any h > 0, we denote by 

(£>)/, = {.x e £>| dist(a;,K 3 \ D) > h}. 

Definition 2.1. Let ft be a bounded domain in M. 3 . We shall say that a portion 
S C <9f2 is of Lipschitz class with constants r > 0, L > 1 if for any point P e S, 
t/iere exzsis a rigid transformation of coordinates under which P = and 

H n Q ro , L = {(a;',a;3) e Q ro ,i|a;3 > V'(^')}- 

w/iere ip is a Lipschitz continuous function in B' rQ such that 

ip(0) = and \\-ip\\c°. 1 (B' ro ) < Lr . 

We say that f2 is of Lipschitz class with constants ro and L if dQ is of Lipschitz 
class with the same constants. 

Remark 2.2. We use the convention of normalizing all norms in such a way that 
all their terms are dimensionally homogeneous. For example: 

IMIc<m(b; ) = U\\L-{B> rQ ) + r ||V^||L~(B; o )- 

Similarly, denoting by D l u the vector which components are the derivatives of 
order i of the function u, 




and so on for boundary and trace norms such as II • II i , , , II • II i , , where 

J 11 M /j"3(an)' 11 M jf~3(aa)' 

is a bounded subset of M 3 whose boundary is smooth enough. 

We will also make use of the following notations for matrices and tensors: for 
any 3x3 matrices A and B we set A : B = Yhj=i and ^ = \{A + A T ). 

By J 3 we denote the 3x3 identity matrix and by Is ym we denote the fourth 
order tensor such that Is ym A = A. 

In the whole paper we are going to consider isotropic elastic materials, hence 
the elasticity tensor C is a fourth order tensor given by 

C(x) = A(x)/ 3 0/3 + 2fi(x)Is ym for a.e. x in ft, (2) 
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where 17 is a bounded domain in M 3, of Lipschitz class, and the real valued 
functions A = X(x) and fi = fj,(x) 6 L°°(f2) are the Lame moduli. We will also 
use Poisson's ratio v(x) — thtt^vt 7 w ■ 

An elasticity tensor C is strongly convex if there is a positive number £0 such 
that, for almost every x in f2, 

C(x)A : A > £o|^| 2 for every 3x3 symmetric matrix A. (3) 

In the isotropic case ([2]), the strong convexity condition takes the form 

n(x) > cto > 0, 2fj,(x) + 3A(a;) > /3q for a.e. x in Q,. (4) 

In these case, the Poisson's ratio has values in an compact subset of (— 1, i). 
More precisely we can estimate 

- 1 + ^£ < v ( x ) < - - ^ for a.e. x in Vl. (5) 
4 2 4 

In the sequel we will make use of the following norm in the linear space of 
isotropic tensors: 

IICHoo = max{||A|| i oc (0) , ||/i||z,oo ( n)} . 

This norm is equivalent to the usual L°° norm for tensors in the space of isotropic 
tensors. 

Our boundary measurements are represented by the Dirichlct to Neumann 
map. As a matter of fact, since we will restrict our measurements to boundary 
data that have support on some subset of the boundary, we will make use of a 
local Dirichlet to Neumann map. 

Definition 2.3 (The Local Dirichlet to Neumann map). Let ft be a bounded 
domain of Lipschitz class and let E be an open portion of dQ. We denote by 
Hcl 2 (Ti) the function space 

Hll 2 {Y<) := {0 e i/ 1/2 (9fi) : supp </> C e} 

and by H C o^ 2 (Yt) the topological dual of Hc/, 2 (S). We denote by < •, • > the 
dual pairing between i?c/ 2 (E) and iJ co 1 ^ 2 (E) based on the £ 2 (E) scalar product. 

1 /2 

Then, given ip G H C o (E), there exists a unique vector valued function u £ 
H 1 ^) weak solution to the Dirichlet problem 

\ div(CS7u) — in 0, 
\ u = if) on dfl. 

We define the local Dirichlet to Neumann linear map Aj? as follows: 

Ag : V G #V a (E) _> (CV U )n| s G H" 1 ^). 
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Note that for £ = dfl we get the usual Dirichlet to Neumann map. For this 
reason we will set A c := Ag n . 

The map A^ can be identified with the bilinear form on Hll 2 {Yj) x ffc/ 2 (£) 

by 

Ag(^, 0) :=< AgV, >= / CVu : V« (7) 

for all ip, 4> £ iJco 2 (S) and where u solves ^ and v is any function such 

that v = on dil. 

(1/2 —1/2 \ 

HZ*(T,),H^''(T,)} defined by 

||T|U = sup{<TV,0> 1^,0 6^(2), ||V|| fl v» (E) = IMIffV=( S) = !> 

for every T e £ (hU 2 (£), F~ 1/2 (£)) . 

Wc finally recall an extension to systems of Alessandrini's identity |Allj . [Is] . 
Let Mi and «2 be the solutions to 

div(C fe Vu fc ) = in Q 
for k = 1, 2 respectively and with C fc , fe = 1,2, satisfying ([3]). Then 

/ (C 1 - C 2 )Vmi : Vu 2 =< (A C i - A C 2)w 2 , «i > (8) 

where Agi, A C 2 denote the Dirichlet to Neumann map corresponding to C , C 2 
respectively. 

2.2 Main assumptions and statement of the main result 

Let A, L, a , (3 , N be given positive numbers such that JV e N, Qo € (0, 1), 
Po G (0, 2) and L > 1. We shall refer to them as the a priori data. Let ro be a 
positive number. 

Our main assumptions are: 

(Al) O C K 3 is an open bounded domain such that 

|0| < Arl 

and we assume that 

where Dj , j = 1, . . . , N are connected and pairwise non overlapping open do- 
mains of Lipschitz class with constants ro, L. 
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We also assume that there exists one region, say D\ such that dDi n d£l 
contains an open flat portion E and that for every j £ {2, . . . , N} there exist 
jl, ■ ■ ■ ,3m S {1, . . . , N} such that 

Dj! = D x , D jM = Dj, 

and, for every k = 2, . . . , M 

oil, n dD jk 

contains a flat portion E^ such that 

E fe c Q, Vfc = 2, . . . , M. 

Furthermore, for k = 1, . . . , M, we assume there exists Pk £ E^ and a rigid 
transformation of coordinates such that Pk = and 

SfcnQ ro / 3 ,£ = {ieQ ro /3,i : a;3 = 0}, 

D jk n Qr /3,£ = {> £ Qr /3X ■ %3 < 0} , 
^Jfc-i n Qr /3.L = {x£ Qr /3,L ■ ^3 > 0}; 

where we set Si := E. 

For simplicity we will call Dj 1 , . . . , Dj M a chain of domains connecting D\ 
to Dj. For any k £ {1, ...M} we will denote by nk the exterior unit vector to 
dD k in P fe . 

(A2) We assume the tensor C to be isotropic piecewise constant of the form 

N 

C = E C ^( a; ) ( 9 ) 

where 

Cj = AjJ 3 (g) J 3 + 2fljlsym 

with constant Lame coefficients Ai and fii that satisfy 

"o < Hj < cxq 1 , X 3 < 2 N + 3A, > fy, j = 1,...,N. (10) 

For j = 1,.,.,N, we denote the Poisson's ratio by v» = — r. Note that 
each Vj satisfies (|5|). 

In the sequel we will introduce a number of constants that we will always 
denote by C. The values of these constants might differ from one line to the 
other. 

Theorem 2.4. Let £1 and E satisfy assumption (Al). Then there exists a 
positive constant C depending on L, A, N,ao, /3o only such that, for any C k , 
k = 1,2 satisfying assumption (A2), we have 

IIC 1 -C 2 ||oo < C||A& - Ag 2 ||*. (11) 
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A better evaluation of constant C is given in Remark 14.71 

In order to prove Theorem 12.41 we will first state it in terms of the forward 
map that maps Lame parameters to the corresponding Dirichlet to Neumann 
map. Then, we will apply to the forward map the following general result: 

Proposition 2.5. Let Mi and M 2 be positive numbers and d E N. Let A and K 

be an open subset and a compact subset ofM. d respectively. Assume that K C A, 

dist (K,R d \A) > Mi , and K C B Ah (0). 
Let B be a Banach space and let F : A^t B be such that: 
(i) F is Frechet differentiable; 

(ii) the Frechet derivative F' : A — > £(R d ,B) is uniformly continuous with a 
modulus of continuity (Ji(-); 

(Hi) i 7 ] K is injective; 

(iv) ( K ) 1 : -F^K) — > K is uniformly continuous with a modulus of continuity 

"2(0; 

(v) F' is injective in K, namely there is a positive number qo such that 

mm \\F'(x)[h]\\ B > 90 ; 

x^K.,\h\ — 1 

then we have 

\\xi - x 2 || R d < C||F(xi) - F(x 2 )\\ B for every x u x 2 G K, 

where C = max{ -gj^-p for 6 1 = ^iam{6 , M 2 } with 6 Q = trf x (f ). 

This proposition holds also in inhnite dimensional spaces. For a proof of 
Proposition 12.51 in finite dimensional space we refer to [BaV, Prop. 5 ]. 

Let us now introduce the forward map corresponding to our problem. In 
order to represent the set of Lame parameters, we will use the following notation: 
let L := (Ai, . . . , Ajv, Mi, ... , /xjv) denote a vector in M. 2N and denote by A the 
open subset of R 2JV defined by 

A := {L e M. 2N : Mi > 0, 2^ + 3Aj > 0, j = 1, . . . , N}. (12) 

For each vector L S A we can define a piecewise constant isotropic elastic tensor 
Cl (as in ^) with Lame parameters Xj and fij for j = 1, . . . , N. 
In this case, ||Cl||oo is equal to the norm in M. 2N given by 

Halloo = . max {max{|Aj|, 

j = l,...,N 

The (nonlinear) forward map can be defined as follows: 
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Definition 2.6. Let £1 and E satisfy assumption (Al). 
Let us define F : A -> £(i?i/ 2 (E), F~ 1/2 (E)) 6y 

We can identify F with F : „4 ->■ B such that = (defined in ©), 

where 2? is the Banach space of bilinear form on h11 2 {Y.) x iJc/ 2 (E) with the 
standard norm. 

1/2 

Let ■0 and </> G -ffco (E) and let be the solution to 

f div(CiVu) = in 17 
1 u = t/» on dtt, 

and solution to 

f div(CiVu) = in 
]_ i) = 4> on 957, 

then 

^ <t>) =< F{L)(ip), (j>>= <C L Vu L : Vv L - 

In the sequel, we will write F and A^ £ instead of F and A|? . 

With the above notation, Theorem 12.41 can be stated as follows: 

Theorem 2.7. Let Q and E satisfy assumption (Al) and let K C A be the 

compact subset 

K:={L€A:a <»j< a^\X < \ 2 Mj + 3A, > /3 , j = 1, . . . , TV}. (13) 

Then, there exists a positive constant C , depending on L, A, N, cto, /3q only such 
that 

W^-L^KCWFiL^-FiL 2 )^ 
for every L 1 , L? in K. 

Notice that Theorem 12.71 means that F is invertible on K and its inverse is 
Lipschitz continuous. 

In Section 2] we will show that the forward map of definition 12.61 satisfies 
all the assumptions of Proposition 12.51 for A and K defined as in (fT2)l and (| 13[) 
respectively and B is the space of bilinear form on iJc/ 2 (E) x £fco 2 (E). Then, 
Theorem 12.71 is a consequence of Proposition 12.51 

3 Preliminary results 

3.1 Further notation and definitions 

In order to prove the main theorem we need to introduce some further notation 
and definitions. 
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Construction of an augmented domain and extension of C. 

First wc extend the domain 51 to a new domain 51o such that 8Qq is of Lipschitz 
class and B ra / C {P 1 ){^Y 1 c ^o, for some suitable constant C > 1 depending only 
on L. We proceed as in |Al-Ro-R-Vel Sect. 6]. Set 



pi = r /C L , where C L = 3 ^ /1 + L ; (14) 



and define, for every x' £ B'ra 

3 



f for < ft, 



P1-2LM for ffc < |^'! < ft, 

for > ^. 

Observe that, for every x' € B' m , \ip + (x')\ < ^ and iVx^^')! < 2i - Next, 

3 

we denote by 

Do = {x = (x\x 3 ) E Qr /3,L I < x 3 < ip + (x')} , 

n = n u A,. 

It is straightforward to verify that 
i) Qq has Lipschitz boundary with constants 3L. 

ii) 

fio 3 Qr /4LC L ,L- 

Let C be an isotropic tensor that satisfies assumption (A2). We still denote 
by C its extension to fio such that C\d A = 2A for every 3x3 matrix A. This 
extended tensor is still an isotropic elasticity tensor of the form 

N 

£ = Y J ^XD 3 (x) (15) 

j=o 

where each Cj for j = 0, . . . , N has Lame parameters satisfying (fpQ|) . 
Construction of a walkway 

Let us fix j G {1, ...N} and let Dj 1 , . . . , Dj M be a chain of domains connecting 
Di to Dj. For the sake of brevity set — Dj k , k = 1, ...M. 

By [Al-R cTR-Vel Prop. 5.5] there exists C' L > 1 depending on L only, such 
that (Dk) h is connected for every k E {1,...M} and every h £ (0,ro/C' L ). 
Denote by 

7 • fr r pi \ . , 

/in = mm < — . — — , — ■ > , lo 

I 6 ' C' L ' 8V1 + 4L 2 J V ' 

where pi is as in (fT4|) . 

Let us introduce the following sets: 
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i) Q(k), k — 1,...,M, is the cylinder centered at Pk such that by a rigid 

transformation of coordinates under which Pk = and belongs to 
the plane {(x',0)}, is given by = Q Pi /al.l- Moreover we denote 

Q(m) = Q(m) n d M -i] 

ii) /C is the interior part of the set (J^o 1 Df, 

hi) /c ?l = uf=o 1 ( D i)hi for ever y h € (°) M; 

iv) 

M-l 

^=K fl ug m u |J g (fe); (17) 
fc=i 

v) 

# = {zeA> I dist(x,dn) > . 

It is straightforward to verify that Kh is connected and of Lipschitz class for 
every h € (0, h ) and that (in a suitable coordinate system) 

^o^VW^S'f)- (18) 

3.2 Existence of singular solutions 

3.2.1 Fundamental solution in the biphase laminate 

In our proof of Lipschitz stability estimates, as in the approach used by Alessan- 
drini and Vessella for the conductivity equation |AV] , a crucial role is played by 
singular solutions for the Lame system. As a matter of fact we are not only in- 
terested in the existence of such singular solutions, but also in their asymptotic 
behavior close to the interfaces. 

In the scalar case, this tool is granted by the existence of Green functions 
and by explicit expressions for solutions in the presence of an interface. 

Unfortunately, the existence of the Green matrix cannot be inferred for el- 
liptic systems with bounded coefficients in dimension 3 or higher. 

Moreover, whereas for the scalar (even complex valued) conductivity equa- 
tion is fairly easy to get a solution of the equation close to a fiat interface by 
using the fundamental solution for the Laplace equation and suitable reflection 
arguments, this seems not to be possible for the Lame system. 

In order to construct singular solutions, we make use of special fundamental 
solutions constructed by Rongved in [R] for isotropic biphase laminates. 

Consider the isotropic tensor 

Cb = C + xr3 + + C _ xr 3 _ 
where C + and C~ are constant isotropic elastic tensors given by 
C+ = A/ 3 ® h + 2fil sym , C" = A7 3 ® I 3 + 2fi% ym , 
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with A and /i and A' and jJ satisfy (p~0|) . Denote Poisson's parameters by v and 



In [R] an explicit formula for a fundamental solution r : {(x, y) : 
K 3 ,a: ^ j/} ->• R 3x3 of 

div (c 6 vr(-,y)) =S y I a , 

is given. Here 8 y is the Dirac distribution concentrated at y. 

This explicit formula is quite involved and some alternative formulations in 
more convenient tensor form have been proposed, for example in [MeRej . 

For the purpose of the present work we need to point out some properties 
of biphase fundamental solution. First of all, it is a fundamental solution, in 
the sense that T(x,y), is continuous in {(x,y) £ I 3 x I 3 : i / i/} and T(x, •) 



for all x G R , and, for every vector valued function 



is locally integrable in 
4> e C5°(M 3 ), we have 



C b VT(;y):V4> = 4>(y). 

Furthermore, for every x,y GM 3 , i^y, we have 

T(x,y) = T{y,x) T , 
C 



(19) 



|r(z,y)|< 



f - y\ 



and, for any r > 0, 



l|vr(-,y)|| L 2 (K3 \ Brto)) 



< 



c 



(rrg)V2 



(20) 



where C depends on ao , /?o only. 

We will also need to use the explicit representation of the some components 
of biphase fundamental solution. In particular, we will make use of explicit 
expression of the third column of T. For x — (x\, x%, X3) with X3 > and 
y = (0, 0, c) with c > 0, from [R] we have, 



T(a;,y) • e 3 
where 



4(1 - v) 




1 











bZ ] + 






V B I3 / 





1 



1 



1 



4(1-1/: 

~3-4i> 2c(x 3 + c) 
a 

[ i?2 

3(3 - Au) 



(21) 



*2 



Rn 



with 



47T/X 



7 log(i?2 + %3 + c) 



/X- /X 



/i + (3 - 4i>)/x' ' 



(22) 
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F 2 (fj,,n',u, v') :-- 



4(l-i/)/i 



(l-2v)(3- 4i/) -2^// 



(23) 



M' + (3 - 4i/)i/ 

J2i = (s? + + (x 3 - c) 2 ) 1/2 , and i? 2 = (x\ +\ +(x 3 + c) 2 ) 1 ' 2 . 

It is worth noticing that, for C + = C _ , the matrix F coincides with the free 
space fundamental matrix for constant isotropic elasticity tensor. 



3.2.2 Singular solutions 

With the aid of the biphase fundamental solution L, let us construct singular 
solutions with singularities in some subset of the domain. Since T is defined 
only in the case of a flat interface, we will need to keep away from interfaces 
that are not flat. 

For this reason we need to set the following notation: let F be the union of 
the flat parts of uf =0 dDj. By fiat parts we intend that they can be represented 
as the graphs of a constant function in at least a ball of radius ^ (as T, k in 
assumption (Al)). Let V = UjL dDj \ T . The set T> contains the non flat 
parts of the interfaces. 

Let C = Y^j=o ^jXDj with tensors Cj satisfying (A2) for all j. Let y <G ilo\D 
and let r = min(ro/4, dist(y, D U dCla)). Then, in the sphere B r (y) cither C 
is constant, C = Cj or, by a suitable choice of the coordinate system, C = 
Cj + (Cj+i — Cj)x{a;3>a} f° r some j = 0, 1, • • • ,N and some a with \a\ < r. Let 

Cj if C = Cj in B r (y), 

Cj + (Cj+i - Cj)x{ X3 >a} otherwise, 

and consider the biphase fundamental solution to 

di V (C y VT(-,y)) = S y I 3 in M 3 . 

Proposition 3.1. Let f2o and C satisfy assumptions (Al) and (A2). Then, 
for y G £Iq\D, there exists a unique function G(-,y), continuous in Q\{y} such 
that 

f CVG(-,y) : = <j>(y), for every e C™(n a ). 
Jn 

and such that 

G(-,y)=0 ondn a . 
Furthermore, if dist{y,VU dClo)) > Tq/c\ for some c\ > 1, 

\\G(;y)-T(.,y)\\ m{no) <Cr^ (24) 

\\G(;y)\\m(a \B A y)) < CVo) -1 ' 2 , (25) 
where C depends ag, @q, A, L and, increasingly, on c\. Furthermore, 

G(x, y) = G(y, xf for every x, y e O \X>. (26) 
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Proof. Let us set 



G(x,y) 



T(x,y) + w(x,y), 



where w is solution to 



div (cv x w{-,y)) 



div((C-C y )Vr(-,y)) in Sl 
-r(-,y) on do . 



(27) 



Since C — C y = in B r (y) and r(-,y) is smooth and bounded in il \B r (y), then 



/ = (C-C„)Vr(-,y) € L 2 (n ), div/ 6 H-^Qo) and T(;y)\ aQo e H l l 2 {d£l ), 



hence this problem has a unique solution w € 7? 1 (r2o). 

Furthermore, if dist(y,2?U 9f2o)) > r o/ c i for some c\ > 1, 

\\w(-,y)\\m(Q ) < C'(||/(- 1 j/)|| H -i(n ) + ||r(-,y)||Hi/2 (9 n )) < C^o" 1 

where C depends on chq, /?o ^4 and ci. This proves Estimate (j2"5) follows 
from ([Ml) and from (12TJ1) . 

The symmetry of function G follows by standard arguments of potential 

theory [E]. □ 



3.2.3 Two useful properties of the biphase fundamental solution. 

Let Cb, Cb be given by 



where C + , C~ and C are constant and strongly convex isotropic tensors whose 
Lame coefficients satisfy ((4]) 

Let Tc b and be the biphase fundamental solutions relative to operators 



Proposition 3.2. For every I, m € R 3 and every y,z G R+, y 7^ z we have 



Let £ > such that B eo (y), B £o (z) C R 3 and B ea (y) n S £o (z) = 0. Since 
Cb = Cfc in R+ , we get trivially 



C b = C+x R 3 +C- X r3, C & = C+ 



Xu 3 , + C x R 3_ j 





Proof. Let us fix 2, m € R 3 and y,z d R+ , y ^ z and denote by 



w = r Ci ,(-,y)/ , w = F Ctj (-,z)m. 




15 



Hence 



/ (C b - C b ) Wu:Wv= [ (C b - C b ) Vw : Vu, for every e £ (0, e ). 

Jr^ Jl3\(B E (i,)uB e (z)) 

(28) 

Integration by parts in R 3 \ {B £ {y) U S £ (^)) yields, for every e £ (0,eo), 

/ C b Vu : Vu = - / (c b W) n-v- (c b Vu) n ■ v, (29) 

jR3\(S e (a)US E (2)) ^8B e (j/) V 1 JdB c (z) V 7 

where n is the outward unit normal to d (B e (y) U B e (z)). Now by definition of 
fundamental solution we have 

lim / (CbVu ) n ■ v = v(y) ■ I and lim / ( C/,Vu ) n ■ v = 0. (30) 
By O and §U§ we have 



In a similar way, we get 



lim / C h Vu : Vu = -v(y) ■ I. (31) 

R3\(B c (y)UB e (z)) 



lim / C fe Vu : Vu = -u(z) ■ m. (32) 

jR3\(B c (y)UB e (z)) 



By J28J), PU, and ((32]) we have 

Cb — Cb) Vu : Vu = u(z) ■ m — v(y) ■ I — Tc{z, y) I ■ m — r^(y, z)m ■ I. 

By (HU) the thesis follows. □ 



Proposition 3.3. Let h, k be real numbers and let H be the fourth order tensor 

W(x) = (hl 3 (g> I 3 +2kl sym ) X R3jx). 
For every I, m £ R 3 and every y, z £ R+, y ^ z, we have 

I HVr Cb (-,y)/ : VT Cb (-,z)m = ( ^-T Cb+m (y, z)m ■ I 
Jr 3 _ V "* 

Proof. Let us fix l,m £ R 3 and y,z £ R:J_, y ^ z. Let to be a positive number 
such that for every t £ (—to, to) the tensor Cb + tW is strongly convex. 
Since H(x) — for every x £ M.+ we have trivially 

4>^ m) (y,z) := [ MVT Cb (;y)l: VT Cb (;z)m= /nVL Cb (•, y) I : Vr Ct (•, z) m. 
Jm 3 _ Jr3 
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Hence, for every t G (— to, to) \ {0}, 

^ m \y,z) = -( ((C b +m)-C b )WT Cb (;y)l:WT Cb (;z)m= (33) 
= \ ((C b +tM)-C b )VT Cb+m (-,y)l : VT Cb (;z)m- 
- f m(vr Cb+m (;y)l-VT Cb (;y)l) :Vr Cb (-,z)m. 

By (|33|) and by Proposition 13.21 we get 

4> {l ' m) (y,z) = j((T Ci+m (y,z))-r Ch (y,z))m-l (34) 
- f e(vr Ct+tH (-,y)/-vr Ci (-,y)z) : Vr Cb (-,z)m. 

Now, by straightforward calculation on the biphase fundamental solution given 
in ^, we have, for every x e Ei, 

lim (vr Cb+tH (x, y) I - Vr Cb {x, y) l) - 0, (35) 

and, for every t € (— to 5 ^o) 5 

H(vr C6+ffl (x,2/)Z- Vr Ci ,(a;,2/)z) : Vr Cb (a;, z) m| < Cmin{|a;|- 4 , 1}, (36) 

where C depends on ao, /3q, y, z, I and m only. By p5|). (|36| . and applying the 
dominated convergence theorem, we get 

lim f H (vr Cb+tH (., y) Z - Vr Cb (., y) l) : Vr Cb (., z) m = 0. (37) 



Finally, by (|3Tj) and (J34J) the thesis follows. 

□ 



3.3 Some estimates for solutions to the Lame system 

In this section we collect some properties concerning solutions to the linearized 
elasticity system with piecewise constant elasticity tensor that will be crucial to 
proe our main result. 

First we state a regularity result for solutions of elliptic systems in composite 
materials from [LNj and |CKVCj . Afterwords we use this result in order to 
obtain Proposition 13. 5[ which is then used to prove a quantitative estimate of 
unique continuation for solutions of systems satisfying assumptions (Al) and 
(A2). 
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Proposition 3.4. Let C + and C~ be two isotropic elasticity tensors with con- 
stant Lame coeffcients X, fj, and A', fjf respectively, satisfying assumption (A2), 
let R > and C5 = C + x R 3 + C~Xr 3 ■ Let v G H 1 (B^) be a weak solution to 



dw(CbVu) — in Br. 

Then, for every multiindex j3' , D@ v € C°(Br) and v <E C°°(B R ). Moreover for 
any S > and k > 



iZ*||B*«|| ioo(gT 5 < C7JT 3 / 2 



|w| 2 da; 



1/2 



(38) 



where C = C(<5, fe, ao, A))- 

Note that, by u € C°° (s]f) 
and ui , has a C°° extension to Bl 

\ B + R 

Observe that, in particular, we have 



\v\\ Lx(S(i _ S)r) < CR-V* 



we intend that v\ _ has a C°° extension to B R 



\v\ 2 dx 



1I/2 



(39) 



Beside this regularity result, the principal ingredients for proving estimates 
of unique continuation are the three sphere inequality, some stability estimates 
for the Cauchy problem and a smallness propagation estimate in a cone. All 
these results holds for constant elasticity tensors and are precisely described in 
the Appendix. 

The flatness assumption on interfaces allows us to get a better estimate 
for the Lipschitz constant. The reason is the fact that solutions to the system 
with piecewise constant elasticity tensor have analytic extension beyond the flat 
interface. This property is stated in the following Proposition. 

Here we use the notation of section l3.2.1l and assumptions (Al) and (A2). 
Proposition 3.5. Let v € i?/ oc (/C) be a solution to 

dw (cVv) =0 in JC. (40) 

Let us fix k £ {0, . . . , M}. Then there exist two positive constants C\ and C , 
depending only on oto, /?o cmd L, such that Vm k can be extended by a function 
v in the set Dk U S fe l 1; where 

eftU D k+1 U S fc : dist (x, B ro/6 (P k ) n Sfe) < (41) 



and 
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Proof of Proposition [ 
It is not restrictive to assume that Pk = 0, belongs to the plane {x% = 0} 
and B+ g , 3 C Dk, ^~/ 3 C -Dfc+i- By Proposition 13.41 we know that v\ Dk G 
C 00 (Di I U Efe) and Wm fc+1 € C°°(-Dfc + i U S^). Since C is a constant tensor in 
each domain, for every j3' £ (NU {0}) 2 , ff/W is a solution to (|40|) . 



Denote by R — ^ and w + := For any xq in S^nS^ , -Df, i> G C° (Bfl(xo)) 

and v \b r M £ C °°(B + r(xo))- 

Let us recall the following Caccioppoli inequality ([BBFM, pag.20]): let u 
be a solution to (|40|) then, for every xo € fl Br 

I \Vu\ 2 < ° [ \u\ 2 (43) 

JB P2 (x ) (P2~Pl) Jb p1 (x ) 

for < P2 < Pi < R where C depends on ao and /?o only. 
Denote by 

<j)(x') = v + (x',0) axidip(x') = ^±(x',0). (44) 

We now prove that <p and V are analytic in n P>r and we estimate the 
derivatives -Df, </> and D^i/) from above on i^/^zo) for every ft' € (NU {0}) 2 
and for every io € fl .Br. 

Starting from (|43|) and using the same iterative procedure followed to prove 
inequality (39) in jBFj we obtain 

2\ N o 

d:,v\'<\c\^\ / \v\ 



/ in/3' 12 ^ / ^ ( 2N \ 1 / , ,2 



|^|=JV J - B fi( ;r o) \ \ 11 / / JBr(x ) 

for every iVo G N, where C depends on ao and /3o only. 
On the other side, by Proposition 13.41 we have 

\\ d I< v \\l°°(b r {x )) < C\\ D %'V\\l*(Br(x ))- 
4 2 

Hence, proceeding as in jBFl (41)] and applying Proposition ^. 41 to D^,v, we get 

\P'\ 

(45) 



f C\ 1 1 

\\D^,(/>h^{B'^ )) ^ I 3 ' 1 ( fl) ll«IU*(-B*(»o))i 



and 

/ r \ 1/8' l+i 

Pf^lU-^ao)) < /3'! (^-j IMU W * o)) . (46) 

where C depends on ao and f3 only. By (1431) and (|4"6) and by the Cauchy- 
Kowalevski theorem we have that the solution v to the Cauchy problem 

div (cVv) = 



v{-,0) = <j) onS t nB R , 
dx 3 1 
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is analytic in the neighborhood H^ij of Ej-nB^. Therefore, taking into account 
(|4"4"]) . u is the analytic extension of w+ in and estimate (j4"2"j) follows. □ 



Finally we state a quantitative estimate of unique continuation. 

Proposition 3.6. Let E\, E\ and h be positive numbers, h < ho, where ho is 
defined in (|16[) . Let v £ Hl oc (JC) be a solution to 

div (cVv^j = in /C, 

such that 

M\ L °°{K ) < £i, 

and 

W ,)|< £l /^Mj 2 forevery xe JC h/2 . (47) 

Then 

it,(*)i<cr(2)%r^w-^, 

wftere re (0, |y), a: = Pm + mi;, 



8 

r 



Jo, 

and m, C , 5, 9, < 9 < I, depend on A, L, ao, /3q and N . 



The proof of the above Proposition is given in the Appendix. 



4 Proof of the main result 

This section contains the proof of the main result that consists in showing that 
the forward map introduced in definition 12.61 satisfies all the assumptions of 
Proposition 12.51 



4.1 Differentiability of F 
Proposition 4.1. The map 

F:A^C{HU 2 ^),H- l l 2 {^)) 
defined in \2. 6\) is Frechet differentiable in A and 

<F / (L){H]iP,0>= [ HVu c :Vwc (48) 

where H = Cjj. 

Moreover, F' : A -> C (R 2 ^,£(#i/ 2 (£),iI^ 1/2 (£))) is Lipschitz continu- 
ous with Lipschitz constant Cf> depending on A, L, ao, $o only. 
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Proof. Fix Lei and let H £ R 2N such that \\HWa, is sufficiently small. By 
|J8]) we have 

< {F(L + H)- F(L)) i=A>= [ BVul+h : Vv L . 

Jn 

Hence, by setting, 

n :=< (F(L + H)- F(L)) ijj, <f> > - [ HVul : V«l = / UV(ul+h-ul) ■ V«l, 

Jn Jn 

we have 

\V\ < Crl\\H\U\V(uL + H-UL)\\ L Hn)H\\ H ij2^y (49) 

where C depends on A,L,ao, /3o only. Let us estimate ||V(ul_|__/j — itj,)||i,2/Q). 
For, observe that w := ul+h — ul is solution to 



(50) 



div(CiVw) = div(HVw_L +i£ ) in SI, 
w — on 9f2, 

hence 

/ ClVw : Vw = / HVul+h : Vi«. 
By Lax-Milgram Theorem and Korn inequality we get 

l|Vw|| L 2 (f2) < C\\H\\oo\\VuL+H\\mn) < Ctq ^HWooWWhI/z^v ( 51 ) 

where C depends on A, L, ao, /3o only. 
By inserting (j5~Tj) into (|49p we get 

hi < CrollSllLllVi^^l^ll^^), (52) 

where C on depends on A, L, ao, /?o on ly that yields P5|) . 

Let us now prove the Lipschitz continuity of F' . Let L} , L? £ A and set 

w := < (F'(L 2 ) - F'iL 1 )) [ffftM >= [ HVw^ : - [ HVu^i : Vv^i 

Jn Jn 

= / H(VMi2 - Vui,i) : Vt)i2 + / HVw£i : (Vu^a - V^i). 
By reasoning as we did to derive (f52"j) we obtain 

|w| < CFTolliilloolli 2 -^IUIIV'IIjjV^sjII^IIjjI/^s)) 

where Cpi depends only on A, L, ao, and /3q. □ 
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4.2 Inject ivity of F\ K and uniform continuity of (F K ) 1 

In the present subsection we will prove Theorem 14.21 whose statement is given 
below. 
Let 




where 5 £ (0, 1) is as in Proposition 13.61 The function a is strictly increasing, 
concave, and lim t _j. &(t) = 0. We have 

Theorem 4.2. For every L 1 ,L 2 £ K the following inequality holds true 

WL 1 - L 2 \\oo < C^iWFiL 1 ) - F(L 2 )\U) (54) 

where o~ N (-) is the composition of the function &(■) defined in f53)) with itself N 
times and C* is a constant depending on A, L, olq, (3q, N only. 

Remark 4.3. Observe that Theorem 14.21 provides the injectivity of F|k and an 
estimate of the modulus of continuity of (-F]k) _1 - 

In order to prove Theorem l4.2l we need to prove first some preliminary results. 
Let j £ {1, . . . , N} be such that 

and let Dj 1 , . . . , Dj M be a chain of domains connecting D\ to Dj. For the sake 
of brevity set Dk = Dj k . Consider K, K and JCh defined as in Section 3. Let 
W fe = IntfU^LoZ^), U k = n \W k , for k = 1, . . . , M - 1. The tensors and 
C L 2 are extended as in (fT5|) in all of CIq. To simplify the notation we will set 

C := C^i and C := C £ 2. Finally let JCk — fCh, n 144 and for y, z £ fCk define the 
matrix- valued function 

S k (y,z):= [ (C-C)(-)VG(-,y) : VG(,z), 

whose entries are given by 

Sj*' q) (y,z):= [ (C-C)(-)VGW(-,j,) : V&"\-,z) p,g=l,2,3 
Ju k 

and where G^ (•, y) and G^ (•, z) denote respectively the p-th column and the 
g-th columnn of the singular solutions of Proposition 13.11 corresponding to the 
tensors C and C respectively. From (f25j) we have that 

\s¥< q) (y,z)\<C(d{y)d(z))-V 2 ,y,z£jC k , 

where the constant C depends on the a priori parameters only and d{y) = 
d(y,Uk), d(z) = d{z,Uk)- For any fixed q — 1, 2, 3 let us denote by S^' q '(-,z) the 
vector valued function whose elements are sjF' q \y, z), p — 1,2,3; analogously 

wc define S^\y,-), for any fixed p = 1,2,3. 
First we prove 
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Proposition 4.4. For all y,z € /Cfc we /iai>e 
belong to Hh c (JCk) and for any fixed q € {1, 2, 3} 



div(CX7S { k ' q) (- 7 z)) = miC fe , 
and /or any /lied p £ {1, 2, 3} 

dw(CV5^ , '' ) (.,«))=0 in^ fc . 



(55) 



(56) 



Proof. For seek of simplicity, in the proof we will omit the index k. Let us fix 
q G {1,2,3} and let us first show that the vector valued function S^' ,q \-,z) g 
Hi oc (K) for fixed z £ K.. Let cj) & Gg° (B r (y )) where yo G 1C and B r (y ) C /C. 
Consider, for fixed p,q £ {1, 2, 3} 



df/. 



SM{y,z)djtt>(y)dy = 

(C-C)(x)V x G<*>(x,y) : V B G (9) (a;, z)d j( f } (y)dx 

! 

Observe now that by (f2"5")l and by the fact that B r (yo) C K, we have 
/ f \\/ x G {p) {x,y)\ 2 dxdy < +oo, / f \\J x G {q) (x, z)\ 2 dxdy < +oo. 

JB r (ya)Ju JB r { V0 )Ju 

Hence, by Schwarz inequality, we have that, for fixed z G /C, 

(C-C)(z)V*GW(z,y) : VsG^foz^fo) e i X (W x S r (y )), 
so that we can interchange the order of integration and get 

S^(y,z)d J <j ) (y)dy = 

(C-C)(x)W x ( G^(x,y)d j( p(y)dy:V x G^(x,z) ,/./■ 



JB T {y ) 

and using the symmetry of G almost everywhere in U, ()26|) . we get 



-(vo) 



G^{x,y)d^{y)dy 



-(vo) 



G ip) {y,x)djcj){y)dy 



d yj G ( P\y,xWy)dy, 



B r (y ) 



so that 



S^"\y,z)8 j <f>(y)dy = 
(C-C)(.x)V 2 



d yj GW(y,x)<P(y)dy: V x G^(x,z) 



B r (y ) 



dx. 
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Now recalling that G(y,x) = w(y,x) + F(y,x), by the properties of T and by 
the boundary value problem satisfied by d Xj w(y,x) for any j — 1,2,3 it is 
straightforward to see that V y d Xj w(y, x) G L 2 (U x B r (y )) and hence 

d Vj S^(y,z) = f (C-C)(x)V x d n G^(y,x):V x &"\x,z)dx. 
Ju 

Now, arguing as in the first part of the proof and considering now a vector- 
valued test function $ G (B r (yo)) and by (|2l)|) we have 



C(y)V v S^(y,z);V y $(y)dy- 



Br(yo) 

(C-C)(as)V» 



s r ( ao ) 



C(y)V y G(y, x) : V y <S>(y)dy 



V x G^(x,z)dx 



and since 



/ C(y) V B G<*> (y, x) : V y $(y)dy = 0, a.e. in W 



for all p = 1, 2, 3 we finally have 



C(2,)V,S^(y,z) : V y <%)dy = 0, 

s r (yo) 

for all $ G C§° (B r (yo)) and since yo is arbitrary ([55j) follows. Analogously we 
get ([56]). □ 



Proposition 4.5. If for a positive Eq and for some k G {1, • • • , M — 1} 

|5 fc (y, z)| < eo^ 1 /° r ewe n/ (y^ 2 ) G -ftTo x iio (57) 

(58) 



|S*(«r,^)|<<>oM- ^ 

V r / V r J \Gi+eo 

where y r = P k+ i + m k+ i,z f = P k +i + fn k+1 , P k+1 G Sfc+i, f,r G (0,r o /C), 

t,. = ^ f;^) > T f = $ (ro) an< ^ n,G,C\,8,6 G (0,1) depend on A, £,ao, /?o 
on??/. 

Proof. Fix z G -Kq and consider the function w(y) := 5£ '^(y, -z), for fixed By 
Proposition 14.41 we know that v is solution to 

div(CVu(0) =0 in t k . 
Moreover, from Proposition 13. 11 we get 

Hy)\ < Giro 1 , yG^ fc , 
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where C\ depends on A,L,oio,/3q only and from ([57)1 we have 

\v(y)\ < r~ 1 e„, yeK Q . 
Then, applying Proposition 13.61 for e\ — Tq 1 £o and E\ — Cir^ 1 , we have 

§" k T r 



\v(y r )\=\S^(y r ,z)\<Cr^{^) 
Now let us consider, for fixed p, 



•£o 



which is solution to 

div(CVu(-)) - in K k 

and which satisfies 

, z€K . 



W«)ISf*c(2) (_ 

By Proposition 13. II we have 



1/2 /jf^\ -1/2 

, z £ K. k , 



where C depends on A, L, a>o, /?o only. Hence, applying again Proposition 13. 61 to 

v(z), we get 

. rail.. «J 



W*)|<Of'(2)""(2 



,Ci+£ 

which proves (|58|) . □ 

Proof of Theorem \4--2\ 
Observe first that \\F(lS) - F (L 2 )||* = || A c - A e ||*. Denote by 

e:=\\F(L l )-F(L 2 )\\*. 

Then from identity ([5]), we derive that for every y, z € Kq and for |Z| = |m| = 1. 



C)(x)VG(x,y)l : VG{x,z)mdx 



< Cr^e, (59) 



where C depends on ao, (3q, A, L only. 
Let 



4 := max {|Aj - Aj-|, \nj - 

0<j<k 
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where k € {0, 1, • • • , M}. Recalling that by construction 

we have that 5q = 0. In order to obtain (|54p we use a recursive argument. 
More precisely, we prove that for a suitable increasing sequence {wfc(e)}o<fc<M 
satisfying e < u> k (s) for every k = 0, . . . , M we have 

4 < ^fc(e) 4+i < w fe+ i(e), for every fc = 0, . . . , M - 1. 

Without loss of generality we can choose ujq(e) = e. Suppose now that for some 
k e {1, • • • , M - 1} we have 

5 k <u) k (e). (60) 

Consider 

S k (y,z):= f (C-C)(-)VG(-,y) : VG(-,z) 
and fix z € l£o- From Proposition 13. II and from (f5T)| we get, for y,z 6 7<o 

where C depends on A, L, ao,A) only. By (|55|) and choosing f — cr with c e 
[1/4, 1/2] we easily get that there are constants Co, <5 € (0, 1) and depending 
only on A, L, ao, /3q and, increasingly, on M, such that for fixed l,mel 3 such 
that III = \m\ = 1, 



-l / r o\ 9/2 / 
\S k (y r ,Zr)m ■ l\ < Cr y—j <;luik{e), 



(61) 



where 



^(M)=( — ) (62) 



Let us choose I = m = e 3 and split 

Sk(.y r ,Zf)e 3 ■ e 3 = I x + I 2 , (63) 



where 



and 



h= [ (C-C){x)VG(x,y r )e 3 : ^G{ Xl z f )e 3 dx (64) 



J 2 = / (C- C)(x)VG(a:,j/ r .)e3 : VG(x,z f )e 3 dx 

Ju k+1 \(B ri nD k+1 ) 

and where r*i = jjaJZ as m 03}- Then, from Proposition 13. 11 we derive 

immediately that 

\h\ < -. (65) 
ro 
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By (f6"U| we have that 

|Afe - Aft | < Wfc(e), |p fe - Mfel < w fc (e). (66) 
and hence, using (fTU|) . 

|vfc-wfc| <CW(e), (67) 

where C depends on ao only. Estimates (j6"6")l and (|67|) together with (pM|) and 
(EH give 



|/il> 



(C£ +1 - c^' +1 )(-)vr fe+1 (., 2/r ) e 3 : vr fe+1 (-,* f ) 



a 3 



b t , nD, 



C Cu k {e) 



/r r 



(68) 

where r^ +1 and T k+ i are the biphase fundamental solutions introduced in Sec- 
tion 13.21 where the elastic phase correspond to elastic tensors C^ +1 and C^ +1 
given by 



c k b +1 
c k b +1 



-fe+i; 



3 + Cfc + lX R 3 



up to a rigid transformation that maps E^+i into £3 = 0. Furthermore by (|61[) . 
(l63l and (l65t we obtain 



(69) 



where C depends on A, L, a , ft only. Hence, by (jBT)1) and (|6"8")l and by perform- 
ing the change of variables x — rx 1 in the integral at the right hand side of 
we have 



t +1 -Cj +1 )( I ')Vr H1 ( I ',e 3 )e 3 : Vr fe+1 (x',ce 3 )e 3 dar 7 



po/ 7 ' 



1*0 



(70) 



where 



Since we have 



«,( -i-c 



r \7/2 



— ? ^fc(e),— 



''n 



(C* +1 -C^ +1 )(x')Vr fc+1 ( a; ',e3)e 3 : Vf fe+1 (a/, Gemote' 

-\s; 0/r 



< c- 



where C depends on A, L, a , (3 only, by ([70)1 and by applying Proposition 13.2 
we obtain 

( T \ T 

(r fe+ i(e 3 ,ce 3 ) - r fe+ i(e 3 ,ce 3 )) e 3 • e 3 < S — + C— . (71) 
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For seek of simplicity in what follows we will omit the indices k and k + 1 
and write /i = fik, n' — /J-k+i and, in a similar way, define A, A', v, v' . We will 
bar corresponding Lame coefficients for C. 

Using the explicit form (f2"T|) of the biphase fundamental solution we have 



rfe+i(e 3 , ce 3 )e 3 • e 3 - r fc+ i(e 3 , ce 3 )e 3 • e 3 
f 



1 



1 1 



+ 



16tt(1 + c) 

c 

47r(c + l) 3 



47r(l — c) \[i /i / 
a [(3 - Av) 2 -7 + 3-4i/] a[(3 - 4P) 2 - 7 + 3 - 4P] 



M(l - * ) 
a 



p(l-P) 



where 



7 = F 2 0u,//,2/,z/), 7 = F 2 (fl,fl / ,D,D') 

and i*i and F 2 have been defined in (|2"2"j) and (|2"3")l . 

From ([7T|) . (|B6l and (j6"7| we obtain, by elementary calculation, for every 
ce [1/4,1/2] 

(72) 



|p(c)|<C U, _ + _ +Wfc ( e ) 



ro 



''0 



where 

p(c) := [4(q - a) (3 - 4i/)(l - u) - (07 - 07)] (1 + c) 2 + 4c(a - a) 

and C depends on ao only. Now, if Wfc(e) < 1/e then we choose r — r £ where 

r E = ^|logw fc (e)|"37, 

where C depends on ocq, /3q, A, L and (5 € (0,1) is as in Proposition 14.51 and by 
([72]) we get 

|p(c)|<C|logw fc (e)|-w, (73) 

for every c G [1/4, 1/2], where C depends on a 0l A)j -<4, £ an d 5 € (0, 1). 
Otherwise, if Wfc(e) > 1/e, since p(c) is bounded, we can trivially write 



W)\ < Cecums) 



Estimates (J73) and ([HI) yields 



|p(c)| < Ca(u k (e)), 



(74) 
(75) 



where C depends on ao, Po, A,L only. 
From ([75]) we easily get 



\a-a\< Ca(uj k (e)) , [cry - cry! < Ca{uj k {e)). 



(76) 
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where C depends on A, L, ao and /3q. 

Now, by simple but tedious calculations, from (l76l) and (J5)) we derive 

W - v>\ < Ca(u k (e)), |V - A'| < Ca(oj k (e)), |// - 0\ < Ca(oj k e)) 

Hence we have 

Sk+i < w fe+ i(e) := Ca(u! k (e)). 
Finally, by iteration and recalling that uio(e) = e we get □ 



4.3 Injectivity of F'(L) and estimate from below of F' K 

Proposition 4.6. Let F : A 

in definition \2.(A Let us define 



Proposition 4.6. Let F : A -)• £(#c /2 (E), i7 co 1/2 (£)) be the map introduced 



go := min (L) [H}\L \L£^H£ R 2N , \\H\L = 1} 

we have 

(^r 1 (1/C*) < «>, (77) 
where cr(-) is defined by (|53p and C*, C* > 1, depends on A, L, ao, /3o AT 

Proof. By the definition of go we have that there exist L £ K and i? = 

(fto.l) • ■ • , ^0,Af, feo.l, • ■ • , fco.Jv), HSolloo = !) sucn tliat 

ll-F' ^o)[^o]L = 9o. 
Therefore, by (|4"5|) and Proposition 14. II we have 



l VG(-,y)l: VG(-,z)m 



12 



c 

< — go for every y,z £ IC , (78) 

7*0 



for every /,m unit vectors of R , where C depends on ao,/3o £ and A only, 
Ho = C_h_ , G(-,y) denotes the singular solution defined in Section [3.2. 21 
From now on vector (0, ho i, . • • , ho n, 0, ko i, . . . , ko n) will still be denoted by 
So- 

Let us fix j £ {1, . . . , N} and let Dj 1 , . . . , Dj M be a chain of domains con- 
necting D\ to Dj, where 

max{|/i j|, |fc 0j -|} = IISoIL = 1 

For the sake of brevity set Di — Dj { , i — 1, ...M and order domains and en- 
tries in _ff accordingly. For every i £ {0, . . . , M — 1} denote by IV|_i(-, •) the 
biphase fundamental solution introduced in Section [3 . 21 where the elastic phases 
correspond to the elastic tensor given by 

C 6 +1 O) = QXrs + C l+ iXR3_(a;)- 
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Here Xq — and /xq = 1- Now, for any i £ {0, . . . , M — 1} let us denote by 



% := max {|/l ,p| , \k),p\} ■ 

0<p<i 

In order to obtain (I77|) . we use a recursive argument. More precisely, we prove 
that for a suitable increasing sequence {wi(<7o)}o<i<M satisfying qo < uii(qo) for 
every i = 0, . . . , M we have 

Vi < ^i(<Zo) %+i < w i+i(<7o), for every i = 0, . . . , M - 1, 
from which, taking into account that we can choose u>o(qo) = qo, we will obtain 

For any i £ {0, . . . , M-l}, let = Int(u;. =0 £) i ), U t = Q\Wi, it = K h nWi 



and, for y,z £ JC, let %{y,z) = {t} p ' q) (y, z)X 
function whose elements are given by 



be the matrix valued 



<p,<3<3 



T?' q \y,z) 



floOVGto(-,i,):VG<«>(-,*), p,q= 1,2,3. 



Moreover, for any fixed q = 1,2,3, let us denote by T} ' q \-,z) the vector 

valued function whose elements are T^ p,q ^ (-,z),p = 1,2,3; analogously we define 

T~i P ' \y, ■)) f° r an Y fixed p = 1,2, 3. 

Let us fix m unit vectors of R 3 . By ([75)1 we have, for every y, z £ Kq, 



\%(y,z)l-m\ < 



l VG(-,y)l:VG(-,z)m 



Sinw, 



&G(;y)l:VG(;z)m 



C C 
< — qo h — m 

ro r 



C 



< — (qo+Ui(q )). 



where C depends on ao t ft, i and A only. 

Arguing similarly to the proof of Proposition 14.51 we have that there exist 
Co, C\ such that for every r £ (0,r o /Co) the following inequality holds true 
(recall r\ — j±- where p\ is defined in (fT4|l) 



\Ti{y r ,z f )\ < Cr - 



l fro 



9/2 



s{ui(qo),r/r ), 



(79) 



where is defined as in (|62p . 
Now we have trivially 



f U VG(-,y r )l:VG(-,Zr)m = %{y r ,z f )l ■ m (80) 
JDi +1 nB ri 

U VG{-,y r )l :VG(-,Zf)m. 



Ui\B ri 
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On the other side we have 



H VG(-,y r )Z : VG(-,z?)m 

U>\B ri 

where C depend on A, L, oto, /3q and M only. 

By the above inequality ([50)1 and ([79]l we have 



C 

< ~, 
ro 



l VG(-,y r )l : VG(',z f )m 



D i+1 nB r 



< 



^(i + (^) 9/ %( Wi ( 90 ),,A„)) 



where C depend on A, L, cxq, (3q and M only. 
Denote by H the tensor given by 

W(x) = (h . l+ il 3 <8> I 3 + 2k . l+1 I svm ) xr3_ (x). 
From Proposition l3.1l we have, for every c G (1/4, 1/2) and choosing I = m = e 3 , 



/ HVr i+ i (x, re 3 ) e 3 : Vr i+ i (x, cre 3 ) e 3 dx 
C f/r Q \i/ 2 /r \ 9 / 2 



< 



'•0 



(?) 



(81) 



where C depend on A, L, a$, Aj and M only. 

Now, performing the change of variables x — r£ in the integral on the right 
hand side of (18111 we obtain 



HVr 4+1 (£, e 3 ) e 3 : Vr i+1 (£, ce 3 ) e 3 



n/r 

< c 



< 



-) 1/2 + P) 7/ %^(,„),^ 

r / V r I \ r 



Therefore, for every g £ (0, 1/Ci), we have 



HVr i+1 (-, e 3 ) e 3 : Vr, ;+1 (•, ce 3 ) e 3 



<c(e 1/2 + e- 7/ M^(9o),e)) 



(82) 

where C depend on A,L,ao,/?o an d M only. Now, if cjj(qo) < e 1 /2 then we 
choose = llogw^qo)! 37 , otherwise if uji(q ) > e~ 1 /2 then we estimate 
from above the right hand side of (|82p trivially. Hence, by Proposition 13.31 we 
have 
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dt 



|t=0 



<Ca(ui(qo)), 



(83) 



where u is defined by (|53[) and C depends on A, L, ao, /3o and M only. 
By explicit calculation from (f2~Tj) . denoting by 

A(t) = A i+ i + th , l+ i, n(t) = /itj + i + tJfeo,i+i) 
A(t) 



"(*) = 



2(A(i) +**(*)) ' 



and by 



a(t) = F!(iM,iJ.(t),Vi), 7(4) = F 2 {}ii^{t),v l ,v{t)) 1 
for Fi and F 2 as in ([22]) and (|23"]l. we get 



1 



1 



17 r C l+1 +ffl( e 3' Ce 3) TO '' 



|t=0 



(1 + C) 3 

Therefore, from (|8"3"|) and ([3]) we find easily 



{ [4(1 - i/i)(3 - 4^)"'(0) + (a 7 )'(0)] (1 + c) 2 + 4ca'(0)} 



K(0)| < Ca(uJi(qo)), 
\{<*y)'(P)\ <Ca(u t (q )). 



(84) 



The first condition of (IMl) gives 

4(1 - i/j)^ 



(/ii + (3 - 4^)^ i+ i) 2 
hence, by recalling ([5]) and (fTOj) . we have 



< Ca(u)i(q )) 



\ko,i+i\ < Ca(u)i(qo)). 
Taking into account (pi>)) . the second equation of (IMl) implies 



(85) 



8(1 - Vi)ihni\ +X 



2(A l+ i + ^+1)2(^+1 + (3 - 4^+i)/Xj) 2 



1 cm+i 



< Ccr(wi((7o)), 



hence, again by © and (TTU1) . 

< Ca{uJi{qo)), 

where C depends on A, L, ao, /3o and M only. 
Therefore 

?7i + l < Ui+i(qo) := Ca(uJi{qo)), 
where C depends on A, L, ao, /3q and M only. 
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Finally, by iteration we get 

1 = VM < Ca M {q Q ) < Ca N (q ) 
and the thesis follows. □ 

Remark 4.7. Observe that the above proposition implies that the Frechet deriva- 
tive F'(L) is injective for every L £ K and hence point (v) of Proposition 12.51 is 
completely proved. Therefore we have 

WL 1 - L 2 |U < CWFiL 1 ) - F(L 2 )\U VL\L 2 e K, 

where C = max{ (oa ffi (Ji) , Mj = m in{^,a }, M 2 = a 2 {-) = 

C*<^('). * = (e^r 1 (1/C*), 5i = ±min{£ ,M 2 } and <5 = 5^7 and we 
recall that C* is the constant that occurs in Theorem 14.21 Cpi is the Lipschitz 
constant of F' introduced in Proposition 14.11 and C* has been introduced in 
Proposition 14.61 



5 Appendix 

For the convenience of the reader, we recall here some quantitative estimates of 
unique continuation. Although such estimates have been proved in the general 
case where the elasticity tensor is of class C 1 ' 1 , here we give the statements in 
the special case we are interested in, namely we assume that 

C = A/ 3 <8 7 3 + 2/dWm, (86) 

where A and /i are real numbers satisfying (j4|). 

The following theorem is an immediate consequence of [Al-Ml Theorem 5.1] 
and standard estimate of smallness propagation jAl-Ro-R-Vel proof of Theorem 
1.10] 

Theorem 5.1 (Three sphere inequality). Let u be a solution to the Lame system 

div ( CVuJ =0 in B fl 

for some positive number f. Then, for every n, T2, 7*3, such that < r% < r 2 < 
?"3 < we have 




where C and 8q, < 9q < 1, only depend on oto, /3q, — and, increasingly, on 
The following theorem has been proved in |M-R] . 



33 



Theorem 5.2 (Stability estimate for the Cauchy problem). Let C be as in (86\) . 

Let u be the solution to the Cauchy problem 

div(cVuj = in B+ /3 , 
u(x',0) = h(x') o°nB' ro/3 , 
= g(x')onB' ro/3 , 

where h E H^(B' ro , 3 ) and g E (B' rg ^ 3 ). We have 



IMIl~(b+ /r ) + r o|| Vu|| ioo(B + j < 



0/3' 



where C and 6\, < 0\ < I, only depend on cxq, (3q. 



In order to state the following result of smallness propagation in a cone 
(Proposition 15. 3|) we introduce some notation. Given z € M 3 , £ e M 3 , \C\ = 1, 
7 G (0, -|), we denote by 

S 3i {x-z)-C 



C (z, (,7) = < 1 € M l ^- ^>cos7 

I \x-z\ 

the open cone having vertex z, axis in the direction £ and width 27 and, for any 
p > 0, we denote by 



C p ( 7 ) =C(0,-e 3 ,7)nQ p , fll , 



where = — ^ — • 

' tan 7 



Let 71,72,73 S (0, f ) be such that 71 < 72 < 73- 
Denote by 



to = 



X 



1 + sin 73 ' 
1 — sin 72 



1 — sin 71 
and, for any t S (0,to]j denote by 

Sfc = X k ~ 1 t, w k = -s k e 3 , k e N, 

(fe) ■ (fe) ■ (fe) • 7 r- M 

rj, = s fe sm7 3 , r\ = s k smj 2 , r\ = s fe sm7i, fc e N. 

Notice that we have 

B (fe+i) (wfc+i) C S (w fe ) C B («,) (lot) C C p (7 3 ), for every k e 
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Let r be a given number such that r £ (0, X^o]- Let ko be the smallest integer 

r 
In 



number such that y k 1 < ■£■ and let 



t := x _(fen_1) r. 

Notice that 

Xto<t< t 

and 

Wk = -re?,, r[ ko) =rsinj 1 . (90) 



Proposition 5.3. Let C be as in (|86|) and let u be a solution to the Lame 
system 

div(cVu\ = in C p (73). 

Assume that 

\u\ 2 <E\ (!)L) 



where e, £" are given positive numbers such that e < E. Then we have 

\u(-re 3 )\ < ^LeK-E 1 -*, (92) 

where 

ho S i\ 

vr = e(- t ) , 

C and 9, 9 e (0, 1) depend on cxq, /?o, 71, 72 and 73 on/?/. 

Proof of Proposition \5.3\ 
Let rf 3 for i = 1,2,3 be as in ((55)) . By Theorem 15. II we have, for every A; G N, 

/ \u\ 2 <c(f \uA ([ \uf) , (93) 

JB^ k) (w k ) \JB^ k) (w k ) J \JB r{k) (w k ) J 

where C and 9, 9 G (0, 1), depend on «0j A), 7ii 72 and 73 only. 
Denote by 

<? k :=E-*[ \u\ 2 . 



„(*0 



Since, by (1591) . B (k) C B (k-i), we have 



o-fe < E~ 2 / |m| 2 , for every k > 2. (94) 

J B Sk-i) ( w k-i) 
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By the second inequality in (pTTj) and by ([§3"]) and we get 

(Tfe < Co k _ x , for every fc > 2. 

Iterating the last inequality and taking into account the first inequality in (|9ip 
we get 

(2 \ " 
y\ , for every fc > 2. 

Now, we choose fc = fco in the above inequality and notice that 

^feo-l 

hence 



f \u\ 2 <Ce 2 ^E 2 ^\ 



'• B J*o)("'*o) 

where C depends on ao, /3o, 71, 72 and 73 only. 
Finally, by (J5D) and by the estimate (|3"9")l , 

^ r,.,,! (-re 8 )j (rsm7i) JB rsinyl (re 3 ) 

where C depends on ao and (3q only and (|92|) follows. □ 



We finally end this appendix by proving our main result on quantitative 
estimate of unique continuation for solutions of Lame system with piecewise 
constant coefficients. In this proof the elasticity tensor C is of the form (|15|) . 



Proof of Proposition \3.£\ 
Denote by C 2 = max |6, ALCl, 4Ci, ^ j and by p 2 = where C L and C\ 
are defined in ([H)) and in Proposition 13 . 51 respectively and h n is defined in (TTH)) . 
Notice that C2 does not depend on ro and that p2 < fg- It is not restrictive to 
assume n\ — 63. Let us denote xq = P\ + nij^pi. We have by fp~8|) 

Moreover, by Proposition 13.51 we have that the function t>| B , (p^^o £i) can be 
extended analytically to a function vq on B' p2 (Pi) x ^—4^-, and 

INIIWb, (Pl wLaiL^ <C(^+ei), (95) 



.£1 ^0 A A 

4 '4ur;; 



where C depends on A, L, a , /?o and TV only. 

Let us construct a chain of spheres of radius p 2 /4 such that the first is 
B P2 /i(xo), all the spheres are externally tangent and the last one is centered at 
Pi + ^-n±. We choose such a chain so that the spheres of radius p 2 concentric 
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with those of the chain are contained in B' po (Pi) x ^— 4^-, if) • The number of 
spheres of the chain is certainly smaller than a constant mi depending on L, ao 
and /?o only. 

By an iterated application of three sphere inequality (|87| with r\ = ^ , r 2 — 
^l 2 - , r 3 = p 2 and by we have 

ll w o|| i2 (B 3p2/4 (P 1+ a m) ) < Cfef ^ + £l ) 1 -C 1 , (96) 

where 0$, < 9q < 1, depends on ao and /?o only and C depends on A,L,aa 
and /3o only. Since B 3f>2 / i (Pi + 4fni) D B P2 (Pi), by (jHU) we have trivially 

lko|| L2(Bp2/4(Pl)) < Cef \E X +£l f- e o\ (97) 

By the above inequality and by Caccioppoli inequality, [BBFMl pag.20], we have 

C 

c 

ro 



iV^oll^^^^^^r^+^i) 1 -^, 08) 



where C depends on A, L, ao and /?o only. 

By (|97|) and (|98|) we get the following trace inequality 

IMU»(B' /16 (F I ))+^o||(CV U o)ni|| ff _ i < Cef '(^r+e!) 1 "^ 1 . (99) 

Now, let us recall the following transmission conditions 

v\d =V\ Dl on Ei , (cVvj m = ((CVv) ^ m on Si. (100) 



Let us denote by £~ /32 (-Pi) = D x n 5^/33 (Pi). By Theorem E2J (|100| . ((99 
(147)) and the Caccioppoli inequality we have 



IN* "^Wn)) + -oil IL-cb-^w, < Ce™? 1 ( El + £l f-^\ 

(101) 

where 8q € (0, 1) and 0\ G (0, 1) depend on ao and /3o only and C depends on 
A, L, ao and /3o only. 

Now, we prove by induction what follows: given 777,2 = — — — - and 

Pk = (64)~ fc+2 , there exist a constant C depending on A, L, ao, Aj an d (increas- 
ingly) on M , such that for every fc £ {1, 2, Af — 1} the following inequality 
holds true 



h\D h || £ - (B - | /32(Pfc)) + r \\W ]Dh || ioo( , 



< 



<C£i ie ° OEi+ei) 1 -^ . (102) 

If fe = 1 then (flQ2|) is proved in (flOTj) . Now, assume that (fTMll holds true for 
k, 1 < fc < M — 2. Let us denote by tj^ the analytic extension of /2 nD k ^° 
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\K-h /2 n Dfe+iJ U (3^ n D k ) where /i , and 3^ are denned by QHfy, 

PT|) and (|4"Tj) respectively. Let us denote Xk = Pk — ff n fe- By the induction 
hypothesis we have trivially 



II^IU^ Pfe+l/64 (, fc ))<^ e ° ( £l+£l )i-^o 1+( 2 , (103) 

where C depends on A, L, o>o, /3q and (increasingly) A; only. 

Let us construct a chain of spheres of radius pfc+i/4 • 64 such that the first 
is Bp h+1 /4,s4,{ x k), all the spheres are externally tangent and the last one is 
centered at Pk+i + ^j7ifc+i- We choose such a chain so that the spheres of 
radius pfc+i/64 concentric with those of the chain are contained in Kh Q n C\ Dk. 
The number of spheres of the chain is certainly smaller than a constant TO2 . 

By an iterated application of three sphere inequality (|87|) with n — ^p-, 
T2 = 3p 4 +1 , r 3 = pk+i and by (I103[) we have 



^L-{B Pk+l/l (P k + l) ) ^ H W fcll^(B 3pfc + l/4 (P fc + 1 + ^i„ fc + 1 )) 



< 



where C depends on A,L, ao,/3o and (increasingly) N only. 

Now, proceeding in exactly the same way followed to prove (|101j) we have 



ll«p fc+1 llL=o(B p - +i/ 3 2 (p fc+1 )) +ro ll Vw l^ +1 ll^(B; fc+i/32 (p t+1 )) 



< 



nk + ln'ni + km 2 fc + l„m 1 4-fc m2 

<Ce\ x (E 1 +e 1 ) 1 -^ 9 ° 

where C depends on A, L, ao,/?o and N only. Therefore, inequality (|102[) holds 
true for every k € {1, 2, M — 1}. 
In particular by (|102j) we have 

ll«l^-JlL~(B; M/32 (P M ))+^||V V | I?M _J| iTO(B - M/32(PA/)) < 

^Cei 1 e ° (£ 1+£l )i-«i e ° . (104) 

Now by a rigid transformation of coordinate under which we have Pm = and 
tim = e 3 w e have Q7 M \ = P' Pl /4L x \®> *t)- ^ n what follows, first we derive an 
error smallness estimate in a ball contained in the cylinder -B pi / 8L x (0, ^) and 
then we apply Proposition [53] in the cone C Pi /8l(73)j where 73 = arctan In 
order to make precise such a step we adopt all the notation introduced before 
such a proposition. Therefore we denote H J3 = 2L, t = §ZT+^n~^' Since 



sin 73 = V4L 1 , +1 we have 



. pi V4L 2 + 1 

to — 



4 1 + V4L 2 + 1 ' 
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Let 71 , 72 be such that 



sin 71 = - sin 73, sin 72 = - sin 73 



and set 



__ 1 - sin 72 _ 4V4i 2 + 1 - 3 
X 1 - sin 71 4\/4L 2 + 1 - 1 ' 
Let r £ (0, \to) be fixed and 

t := x-^-^r 

where kg is the smallest integer number such that x k ~ x < — . Moreover let 
P = Pm + tnM- Proceeding as before, by (|104j) we get 

„M-lnm 1 + (M-l)m 2 M _ 1 m , _|_ ( M - 1) m 2 

where C depends on A, L, cto, /?o and M only. 

Now denote by p — min{pjy//4 • 64,isin7i}. By applying Theorem 15.11 with 
r i = P, r 2 — t sin 71 and r 3 = t sin 73 we get 

NUb^cp)) < ^ ^Cs^{E l+ s 2 ) l - e \ 

where C, 8 2 , < 9 2 < 1 depend on L, ao and /3q only. Now we use Proposition 
15.31 and we get 

\v(x)\ <c(^y~* ef{E l + e 2 f-^, (105) 

where x = Pm + rnja, C depends on A, L, eta, A) and (increasingly) N only 
and 

*7r = 2 (9 ■ (106) 
Finally, denoting by = min {#1 , #2, #3} and m = max {mi 4- 1, 7712 + 1} we get 

and by t > xto the thesis follows. □ 
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